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ABSTRACT 

Using a large set of high resolution numerical simulations incorporating non-equilibrium molecular hydro- 
gen chemistry and a constant source of external radiation, we study gas collapse in previously photo-ionized 
mini-galaxies with virial temperatures less than 10 4 K in the early universe (redshifts z = 10-20). We confirm 
that the m echanis m of positive feedback of ionizing radiation on star formation in mini-galaxies proposed by 
iRicotti. Gnedin. & Shul]lJ200 2) can be efficient despite a significant flux of metagalactic photo-dissociating ra- 
diation. We derive critical fluxes for the Lyman- Werner background radiation sufficient to prevent the collapse 
of gas in mini-galaxies as a function of the virial mass of the halo and redshift. In our model, the formation of 
mini-galaxies in defunct H II regions is most efficient at large redshifts (z > 15) and/or for large local gas over- 
density 6 > 10. We show that non-equilibrium chemistry plays an important dynamical role not only during the 
initial evolutionary phase leading to the gas becoming gravitationally unstable inside the mini-halo, but also at 
the advanced stages of the core collapse, resulting in efficient gas accretion in the core region. We speculate on 
a possible connection between our objects and metal-poor globular clusters and dwarf spheroidal galaxies. 
Subject headings: early universe — galaxies: dwarf — galaxies: formation — methods: TY-body simulations 
— astrochemistry 
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1. INTRODUCTION 

The formation of the first (small) galaxies in hierar- 
chical cold dark matter (CDM) cosmologies has been 
a subj ect of intensiv e resea rc h. Both non-co smological 
(e.g. iTegmark et all Il997t iFuller & Couchmanl 1200(1 
Bromm, Coppi, & Larson 2002 ) and cosm ol ogical 



(e.g. Abel et al. 



lAbel. Brvan. & Normanl [2002; 

Yoshi daet al.l2003l) approaches produced a consistent picture 
of first stars in the universe forming inside mini-halos with 
virial masses < 10 6 M Q at redshifts > 20. As the virial 
temperature of such small halos is below 10 4 K, rendering the 
atomic hydrogen line cooling inefficient, the primary coolant 
responsible for the collapse of the zero metallicity gas is 
molecular hydrogen. 

lHaiman. Abel. & Reesl d2000l) argued that soon after the 
very first stars in the universe are born in rare 3-4 sigma 
peaks, the Lyman- Werner photons produced by these stars 
dissociate molecular hydrogen making further small galaxy 
formation impossible. As a result, the main bulk of stars in 
the early universe should have formed later, when more mas- 
sive galaxies with T w \ r > 10 4 K, capable of cooling through 
atomic hydroge n lines, are formed . Cosmological simula- 
tions of Machacek. Brva n. & Abell J2001I) and Yoshid aetall 
(2003) corroborated the above picture by showing that even 
small flux F ~ 10~ 23 ergs s" 1 cm" 2 Hz" 1 sr" 1 of the Lyman- 
Werner radiation is sufficient to suppress the formation of 
mini-galaxies in the e arly universe. 

Ricotti et al. (2002) demonstrated that to gain better under- 
standing of early star formation in the universe one has to use 
cosmological simulations with full radiative transfer. They 
uncovered an important mechanism of positive feedback of 
ionizing radiation from first stars on subsequent star forma- 
tion in mini-galaxies. In their simulations, cosmological H II 
regions formed around the first mini-galaxies in dense fila- 
ments become sites of intense H2 production after the source 
of the ionizing radiation is turned off. In the defunct cos- 
mological H II regions gas cools and recombines at the same 
time, resulting in enhanced non-equilibrium ionization frac- 



tion, and, as a consequence, inc reased production o f molec- 
ular hydrogen. In the model of Ricotti et al. (2002) the flux 
of the background Lyman- Werner radiation is calculated self- 
consistently by performing full radiative transfer from the 
sources of radiation (star-forming regions) inside the com- 
putational volume. The authors showed that the significant 
flux of H2 photo-dissociating radiation does not prevent mini- 
galaxies being formed in and around cosmological H II re- 
gions. 

1 1 1 Fig. 7 of Ricotti et al. (2002) and the corresponding 
mpeg animation in the electronic edition of the publication 
one can see some details of the above process. In this figure, 
a mini-galaxy is seen to form first stars at z < 17 and then to 
undergo a few consecutive star bursts due to positive feedback 
of the ionizing radiation. As one can see from this figure, the 
mini-galaxy is located inside a relatively dense cosmological 
filament. The temperature inside the H II regions is close to 
the canonical val ue of ~ 10 4 K . 

More recentlv. lO'Shea et alJ ( 120051) presented their simula- 
tions of the formation of a primordial star within a region ion- 
ized by an earlier nearby star. T hese semi-cos mological simu- 
lations confirmed the results of Ricotti et al. (2002) by show- 
ing that there is a positive feedback due to photo-ionization 
on subsequent star formation in mini-galaxies. 

Cosmological simulati ons w ith full radiative transfer, such 
as those of IRicotti et alJ d2002). are very computationally ex- 
pensive, as one has to achieve both high resolution (to resolve 
collapsed gas in mini-galaxies) and large spatial coverage (to 
have the rare sources of ionizing and Lyman- Werner radiation 
included in the computational box and to beat the cosmic vari- 
ance). Instead, in this paper we use non-cosmological high- 
resolution simulations to investigate the mechanism of posi- 
tive feedback of ionizing radiation on star formation in mini- 
galaxies in greater detail. We model both non-equilibrium 
chemistry and the impact of the external photo-dissociating 
radiation on collapse of gas in mini-galaxies inside defunct 
cosmologic al H II regions at r edshifts 20- 10. We confirm 
the result of IRicotti et al.l (120021) that the positive feedback is 



2 



Mashchenko et al. 



TABLE 1 

Reactions 



NN 


Reaction 




Rate reference 


1 


H + <?~ — » 


H + + 2e~ 


1 


2 


H + + e~ - 




1 


7 


H + (T -» 




2 


8 


H + fT - 


» H 2 + tT 


1 


11 


H 2 + H+ - 


^H++H 


1 


12 


H 2 + e~ - 


-> 2H+<r 


1 


13 


H 2 + H- 


> 3H 


1 


14 


H"+e" - 


4 H + 2e" 


1 


15 


H- + H- 


*2H + e- 


1 


16 


H- + H+ - 


^2H 


1 


17 


H~ + H+ - 


-^Hj + e- 


1 


23 


H- + 7- 


H + e- 


3 


27 


H 2 + 7 - 


H* ^2H 


1 



NOTE. — Reaction numeration is from Ab el et al] 119971) . 
References. — ( ll lAbel et all 19971 <2, IHutcMnsfi976l (3) this work. 



efficient even for relatively large fluxes of external Lyman- 
Werner radiation with F > 10~ 21 ergs s" 1 cm" 2 Hz" 1 sr" 1 . 

2. MODEL 

Our goal is to simulate the evolution of a small, fully 
photo-ionized portion of a dense cosmological filament con- 
taining one mini-halo in the early universe (z ^ 10) after 
the source of the ionizing radiation is turned off. The gas 
remains exposed to the external source of radiation with 
the photon energies below 13.6 eV (which includes Ly man- 
Werner photons capable of dissociating the molecular hydro- 
gen). We assume the gas to be optically thin. To run the 
simulation s, we use an adaptive-mesh P 3 M-SPH code Hydra 1 
dCouchman. Thomas. & Pearcel 19951) . We use the scalar ver- 
sion 4.0 of the code, with the m olecular h ydroge n chemistry 
module developed by Fuller & Couchman (2000). 

2.1. Chemistry 

Molecular hydrogen is the primary coolant in the gas of 
primordial composition inside small dark matter (DM) halos 
with virial temperature less than 10 4 K. In the primordial dust- 
free gas molecular hydrogen can be produced through two 
main channels: Hj channel and H" channel. In our analy- 
sis, we consider only the latter channel, as it is known to be 
the dominant one at redshifts z < 200. In the H" channel, the 
electron acts as a catalyst: 

H + e"^H" + 7 
H" + H^H 2 + e" 

At the temperatures considered here (T ^ 10 4 K) helium 
is neutral. As a result, we only have to follow the evolu- 
tion of five different species: H, e", H + , H", and H 2 . We use 
all 13 major chemic al reactions invo lving these species from 
the compilation of A bel et aD (119971) . The reaction s are listed 
in Tab le [0 We keep the reaction numer ation of |Abel et al J 
( 119971) . Most of the reaction rates are from lAbel et alJ d!997). 
with the two excepti ons: (1) for the reaction 7, we use a sim- 
pler expression from Hutchins ( 1976, this is the way this reac- 
tion is implemented in the publicly available Fortran routine 
on Tom Abel's primordial gas chemistry page 2 ): 

1 http : / /coho .mcmaster . ca/ hydra/ hydra_consort ■ html| 

2 http : / /cosmos . ucsd . edu/$\sim$tabel/PGas/| 



h = 1.83 x io _18 r 08779 s _1 , (i) 

where T is the gas temperature in K; (2) for the photo- 
detachment cross-section of H" (reaction 23) we used our 
own (higher q uality) fitting formula to the original data of 
IWishartU T979) tabulated for photon energies 0.76 . . .9.9 eV. 
More specifically, instea d of the fitting form ula a^i = 7.928 x 
IOV-j/o) 1 V 3 cm 2 of lAbel et alJ dl997h . we used the fol- 
lowing one: 

023 =2.241 x 10V-i/o) L V 31 cm 2 . (2) 

Here frequencies are in Hz. Both expressions are valid for 
hv > hvQ = 0.755 eV. The standard an d maxi mum deviations 
of our fit from the data of IWishartl (119791) are 0.007 and 
0.04 dex, respectively, which is much better than for the orig- 
inal fit (0.03 and 0.07 dex, respectively). 

The external source of non-ionizing radiation is considered 
to have a power-law spectrum: 

^2i=^(^)"l0 21 , (3) 

where Fy is the radiation flux in ergs s" 1 cm" 2 sr" 1 Hz" 1 at the 
average Lyman- Werner frequency with hv = 12.87 eV. The 
rate for the reaction 23 can be obtained from 

K23 = 47T / dv S . (4) 

J hv 

Here hvi = 13.6 eV is the hydrogen photo-ionization energy. 
We consider a steep external radiation spectrum with a = — 1. 
The rate equation is derived by numerically solving the inte- 
gral in equation (0} with 023 (^) given by equation Q: 

fc 2 3 = 7.465 x 10" I0 F 2 i s" 1 . (5) 
The rate for the reaction 27 is 

k 21 = 1.382 xlO" 12 /^" 1 (6) 

(lAbel et alll997l) . 

For a species A, we express the fractional abundance /a in 
terms of the proton number density n of all hydrogen species: 
f A = n(A)/n. Here 

n = p(l -y)/m p = n(H) + n(H + ) + 2«(H 2 ). (7) 

(The abundance of H" is much lower than the abundance of 
other hydrogen species and was neglected in the above equa- 
tion.) Here p is the gas density, m p is the mass of a proton, 
and y is the helium mass fraction. As the destruction of the 
ion H~ proceeds much more quickly than its production, the 
equilibrium fractional abundance of H" can be estimated as 

^ k 7 fnfe- ^ 

(kg + h 5 )f H + (ki(, + kn)fw + kufe- + hi/n ' 

Furthermore, fip = f e - due to neutrality of helium and very 
low abundance of H". In addition, the neutral hydrogen 
fractional abundance can be estimated from equation Q as 
fu = 1 — /h+ _ 2/h 2 ■ This leaves us with only two species whose 
fractional abundance has to be derived by solving differential 
equations: H + and H 2 . The corresponding equations are 
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--kife-fu-kifn+fe--, 



(9) 



n dt 

~^T= hfy-M-M (hifa+kufr + knfa + — VlO) 
n at \ n J 

2.2. Cooling 

Molecular hydrogen molecules cool through radiative de- 
excitation of ro-vibrational l evels. We adopted the cooling 
function from Galli & Pallal ©98 ): 



A H2 = (Alxe + A^q) n H2 «H ergs cm 3 s 



(11) 



where 



At TP. — 



1 



9.5 x 10 22 r 3 3 - 76 ^_ (0 .13/7-3)3 



1+0.127^ 1 



3 x l(T 2 V - 51 / r3 +6.7 x 10- 1 V 5 86 / r3 + 
lixirV 11 - 7 ^' 



(12) 



is the lo cal thermodynamic equilibri um cooling function 
given bv iHollenbach & McKed ( fT979l a sum of their eqs. 
[6.37] and [6.38]) and 

log A„^ = - 1 03 .0 + 97 .59 log T - 48 .05(log T) 2 + 

10.80(logr) 3 -0.9032(lo g r) 4 (13) 

is the de nsity ind ependent low-density limit of the H2 cooling 
function dGalli & Pallall998l) . Here T 3 = Tj 10 3 K. 

We also included fi ve other co oling processes (the rates are 
taken from Anninos et alJll997l) : atomic hydrogen line cool- 
ing, 

A m = 7.50x lO-^l + Tg'/y V 1I8348 / r « e « H , (14) 
hydrogen collisional ionization cooling, 

A ion = 2.18x lO- n hn e n H , (15) 
hydrogen recombination cooling, 

A rec = 8.70 x 10- 27 r I/2 r 3 -°' 2 (l + r 6 °' 7 )~V«H+, (16) 
Bremsstrahlung cooling, 



A b rem= 1.43 x 10- 27 r'/ 2 X 

l.l+0.34 e -( 5 - 5 - lo s r > 2 / 3 l n e n H+ , 



(17) 



and Compton cooling, 



A c = 1.017 x IO-^T^CT-TcmbK. (18) 

Here T$ = T/10 5 K, = T/10 6 K, £1 is the rate for the reac- 
tion 1 (see Table[Q, and 7cmb = 2.73(1 +z) K is the temper- 
ature of the cosmic microwave background radiation, where 
Z is th e redsh ift. The units for the cooling functions in equa- 
tions (11411 8i are ergs cm" 3 s" 1 . The total coo ling functi on is 
the sum of the cooling functions in equations (111 1141 1 81 . 



2.3. Generating Initial Gas and DM Distributions 

In our model, we use periodic boundary conditions with 
the box size 5.7 times larger than the virial diameter of the 
DM halo. As the halo is located inside a relatively dense fila- 
ment, we ignore the expansion of the universe. Both DM and 
baryons have zero global angular momentum. 

The DM halo is assumed to have a Navarro-Frenk- White 
HSTFW: [Navarro. Frenk. & White! 199 7) universal density pro- 
file, 



pDM(r) = 



Pdm.o 



r/r s (\ + r/r s f 



where 



Pdm.o '■ 



ln(l+c)- 



1+c 



(19) 



(20) 



Here r s is the scaling radius of the halo. In the ACDM cos- 
mology, the typical concentration c = r vu -/ r s of a DM halo with 
the virial mass m sa < 10 1 1 M fl at the redshift z is 



27 / m vir 
l+z \10 9 M G 



-0.08 



(21) 



( Bull ock" et alJ200j1) . Here r vlr is the virial radius of the halo. 
Zhao et al. (2003a b) showed that in cosmological TY-body 
simulations c is always larger than ~ 3.5. 

Initially, the gas is isothermal with temperature Tq = 10 4 K, 
with fully ionized hydrogen and neutral helium and is in hy- 
drostatic equilibrium inside the periodic computational box 
containing one DM halo. Our choice for the initial tempera- 
ture of the photo-ionized gas was guided by the r esults of the 
cosmological simulations of Ricot ti et"aT1 i|2[002) incorporat- 
ing full radiative transfer and primordial chemistry. (In par- 
ticular, from their Fig. 7 one can see that the temperature of 
the fully photo-ionized H II region surrounding a mini-galaxy 
is very close to 10 4 K, despite the fact that the sources of 
ionizing radiation in this particular model are Pop III stars 
which have a relatively hard spectrum.) The gas in our model 
has a primordial composition, with the helium mass fraction 
y = 0.24 and no metals. We produce this initial configuration 
in a two-step procedure as described next. 

First, the computational box is filled homogeneously with 
N g SPH gas particles with the temperature Tq. We place a DM 
halo, represented by A^dm = N g particles, at the center of the 
box. The halo is truncated at a radius equal to the half-size of 
the box (5 . 7 times larg er than the virial ra dius). The procedure 
outlined in Mashchenko & Sills (2005a) is used to generate 
the initial coordinates and velocity vectors of the DM particles 
(isotropy is assumed). Then we evolve the cube using Hydra 
for a very long time Ti with the DM particles rigidly fixed in 
their original positions (static DM halo), the gas temperature 
fixed at To (isothermal equation of state), and no chemical 
reactions. We fix the evolution time at T\ = 10 Gyrs for all 
our models, which is ~ 20 - 50 times longer than the sound 
crossing time in the box. At the end of the evolution, the 
gas density profile is well converged to its final hydrostatic 
shape down to the smallest resolved radius (containing 32 gas 
particles). 

Second, we run the simulation for an additional 0.5 - 1 Gyr 
with the live DM halo (gas is still isothermal and fully ion- 
ized) to allow the halo to (slightly) adjust itself to the presence 
of gas in the box. We stop the simulation before the central 
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FIG. 1 . — Initial hydrostatic radial gas density profiles for all our models. Models with larger 5 are located higher, (a) Halos at fixed redshift z = 15. Dotted, 
dashed, solid, and dash-dotted lines correspond to the models with logm v j r = 6, 6.5, 7, and 7.5, respectively. For the case of logm v j r = 7.5, only the models 
with 5 = 1 and 3 are shown, (b) Halos with fixed virial mass m v j r = 10 7 Mq. Dotted, solid, and dashed lines correspond to the models with z = 10, 15, and 20, 
respectively. 



DM density profile starts getting flatter due to two-body inter- 
actions between DM particles. 

One of the free parameters in our model is the average bary- 
onic overdensity in the box, 6 = po/ pb, where po is the average 
gas density in the box and 



TABLE 2 
Halo Parameters 



(22) 



is the average baryonic density in the universe. (Here G 
is the gravitational constant.) For the rest of this paper we 
assume a flat universe with f2 m = 0.27, fib = 0.044, H = 
71 km s -1 Mpc -1 , and cr 8 = 0.84 (Sner gel et alJl2003l) . The 
parameter 8 is meant to represent both the gas overdensity 
and the DM overdensity inside the box. Due to the small size 
of the box (~ 5 proper kpc) the DM which is not a part of 
the virialized DM halo is assumed to be distributed almost 
homogeneously and hence to have a negligible impact on the 
gas evolution inside the box. We do not include this DM in 
our simulations. In our models, when we vary the overdensity 
S, we change only the average gas density, keeping the DM 
distribution unchanged. 

Overall, we have four free parameters in our model: red- 
shift z, DM halo virial mass m v h, overdensity S, and exter- 
nal radiation flux Fi\, By running ~ 150 models we ex- 
plore a wide range in the initial parameters, with z= 10 . . . 20, 
m vh -= 10 6 . . .3 x 10 7 M , 5= 1...30, and F 2 \ = 0. . . oo. 

2.4. Numerical Simulations 

We simulated 6 different DM halos (listed in Table|2j. The 
main emphasis was on halos formed at z = 15 (4 models), but 
we also considered one high-redshift (z = 20) and one low- 
redshift (z = 10) case for the fiducial virial mass of 10 7 M Q . 
Each halo was simulated with four different values of over- 
density 6: 1 (not inside a filament), 3 (low density filament), 
10 (filament), and 30 (partially collapsed halo). Each of the 
above 24 models was simulated with 4-8 different values 
of the external Lyman- Werner flux Fi\ (including the case of 



z log A/lvir C Fs Tvir Tvir Box size £ 

Mq kpc kpc K kpc pc 

li) 7 3.50 0.0964 0.337 4890 3^81 6~6~ 

15 6 3.50 0.0587 0.206 800 2.32 4.0 

15 6.5 3.50 0.0862 0.302 1730 3.41 5.9 

15 7 3.50 0.127 0.443 3720 5.00 8.7 

15 7.5 3.50 0.186 0.650 8020 7.34 13 

10 7 3.55 0.182 0.644 2560 7.17 12 

NOTE. — Here r v j r = AtV y 2 ir / (2k) is the virial temperature of the halo, where 
Vy ir = Gmvir/f"vir, k is the Boltzmann constant, and the mean particle mass is 
/i = 0.633m p for primordial gas with fully ionized H and neutral He; e is the 
gravitational softening length. All units are proper (not co-moving). 



zero flux), which resulted in ~ 120 models ran with full chem- 
istry and cooling. The particular values of F21 (listed in Ta- 
ble|3jl were chosen to bracket the minimum value of the flux 
^2i,crit when the external radiation prevents the halo from col- 
lapsing. 

The virial temperature of our halos is below the initial gas 
temperature of 10 4 K for the case of fully ionized hydrogen, 
and would be above 10 4 K, if the gas were neutral, only for our 
most massive halo with m wlI ~3x 10 7 M Q (see Table|2j>. The 
gravitational softening length e was set to be commensurable 
with the average inter-particle distance within the virial extent 

1 /3 

of the DM halo: e = 1.8rf,/N^ M . (Here r/, is the half-mass 
radius of the halo.) All of the above models had a total number 
of particles of 2 x 64 3 ~ 520,000. 

As one can see from Table|2] the conc entration of our halo s 
is close to 3.5 due to the constraint of Zha o et al . (200333). 
The virial radius of these objects is very small: between 200 
and 650 pc. The size of the periodic computational box is 
between 2.3 and 7.3 proper kpc. 

The central hydrogen number density n c in the initial hy- 
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TABLE 3 
Results of the Simulations 



z 


logm v j r 


5 


l°g^21 


n c 


n c/ n out 


logf^Lcrit 




'"vir 1 


m vir. 1 


r vii,2 


m vir,2 




M 




10~ 21 ergs s cm -2 sr~' Hz~' 


cm -3 






Myr 


pc 


M 


pc 


Mq 


20 


7 


1 


-00,-2,-1,0,0.5,1,2,00 


0.041 


28 


0.75 


31 


7 


7.1 x 10 3 


8 


3.8 x 10 s 


20 


7 


3 


—oo,—l,0,l,oo 


0.13 


30 


oo 


19 


8 


9.0 x 10 3 


8 


9.5 x 10 s 


20 


7 


10 


-oo,-l,0,l,oo 


0.52 


37 


oo 


12 


10 


1.5 x 10 4 


8 


1.9 x 10 6 


20 


7 


30 


—00,0,1,2,00 


4.2 


110 


oo 


4 


263 


1.8 x 10 6 


8 


5.4 x 10 6 


15 


6 


1 


-oo,-3,-2,-l 


1.2 x 10~ 3 


1.6 


— oo 










a 


15 


6 


3 


-00,-5,-4,-3.5,-3,-2,-1 


3.5 x 10~ 3 


1.6 


-3.75 


225 


4 


2.8 x 10 3 


5 


1.6 x Id 4 


15 


6 


10 


-00,-4,-3,-2.5,-2,-1,0 


0.012 


1.6 


-2.75 


102 


5 


3.7 x 10 3 


5 


2.6 x 10 s 


15 


6 


30 


-00,-2,-1.5,-1,0 


0.035 


1.6 


-1.75 


60 


5 


4.5 x 10 3 


5 


6.0 x 10 s 


15 


6.5 


1 


-00,-3,-2.5,-2,-1 


2.0 x nr 3 


2.7 


-2.75 


188 


6 


4.5 x 10 3 


7 


4.6 x 10 4 


15 


6.5 


3 


-00,-3,-2.5,-2,-1 


6.0 x 10~ 3 


2.8 


-2.25 


105 


7 


5.5 x 10 3 


7 


3.0 x 10 5 


15 


6.5 


10 


-00,-2,-1.5,-1,0 


0.020 


2.8 


-1.25 


56 


7 


7.0 x 10 3 


7 


7.9 x 10 s 


15 


6.5 


30 


-00,-2,-1,-0.5,0 


0.063 


2.9 


-0.75 


37 


8 


9.1 x 10 3 


7 


1.9 x 10 6 


15 


7 


1 


-00,-3,-2,-1.5,-1 


7.4 x 10~ 3 


11 


-1.25 


92 


10 


9.6 x 10 3 


10 


3.9 x 10 s 


15 


7 


3 


-00,-2,-1,-0.5,0 


0.023 


1 1 


-0.25 


47 


10 


1.1 x 10 4 


10 


9.9 x 10 s 


15 


7 


10 


-oo,-l,0,l,oo 


0.082 


12 


oo 


28 


13 


1.6 x 10 4 


10 


2.5 x 10 6 


15 


7 


30 


—00,-1,0,1,00 


0.33 


17 


oo 


18 


18 


3.1 x 10 4 


10 


6.3 x 10 6 


15 


7.5 


1 


-oo,-l,0,l,oo 


0.22 


370 


oo 


18 


17 


2.4 x 10 4 


15 


1.2 x 10 6 


15 


7.5 


3 


—00,-1,0,1,00 


0.82 


490 


oo 


11 


27 


5.8 x 10 4 


15 


3.0 x 10 6 


15 


7.5 


10 


















b 


15 


7.5 


30 


















G 


10 


7 


1 


-oo,-3,-2,-l 


1.1 x 10- 3 


4.8 


— oo 


447 


12 


1.1 x 10 4 


14 


3.2 x 10 4 


10 


7 


3 


-00,-3,-2.5,-2,-1 


3.3 x 10" 3 


4.8 


-2.25 


201 


13 


1.3 x 10 4 


15 


7.6 x 10 s 


10 


7 


10 


-00,-2,-1.5,-1,0 


0.011 


5.0 


-1.25 


86 


15 


1.8 x 10 4 


15 


2.5 x 10 6 


10 


7 


30 


-00,-2,-1,0,0.5,1,2,00 


0.038 


5.6 


0.5 


52 


19 


2.8 x 10 4 


15 


5.8 X 10 6 



NOTE. — Here logF2i = — oo corresponds to no external radiation, logftl = oo stands for no H2 cooling, n c and n oul are the initial gas density at the center 
of the halo and at the edge of the computational box, respectively, and fil.crii is the critical value for Fi\ (such that models with F21 ;J ^2l,crit will not collapse 
in one local Hubble time). Parameters t co u, r yII j, m v ; r j, r v ; r 2, and m v ; r 2 correspond to models with F21 = 0. Here ? co u is the time when the Jeans criterion is met 
inside the halo, r v j r 1 and m v j r 1 are the radius and the enclosed gas mass of the region with the smallest virial ratio at t = ? co n , and r v j r 2 and »i v j r 2 are the same 
quantities at the end of simulations. 

a The model does not collapse even when Fi\ =0. 

b The model collapsed during the evolution of isothermal gas in a live DM halo (the second step in generating the hydrostatic initial conditions). 

The model collapsed during the evolution of isothermal gas in the static DM halo potential (the first step in generating the hydrostatic initial conditions). 



drostatic state spans a wide range: 0.001 . . .4 cm" 3 (Table[3j. 
This excludes the two most massive models (logm v i, = 7.5) 
with the largest overdensities, <5 = 10 and 30, which collapsed 
durin g the generation of the initial gas and DM distributions 
(see § I2.3> . As one can see from Table[3] the gas density con- 
trast n c /n out between the center of the halo and the edge of 
the box initially ranges from 1.6 (i.e. the gas is barely per- 
turbed by the shallow gravitational potential of the DM halo) 
for our lowest mass halo with m v j r = 10 6 M Q to ~ 500 for the 
model with m wu ~ 3 x 10 7 M Q and 8 = 3. For models with 
larger m Yn , 8 and/or z the ratio n c /n oM increases noticeably 
with 8, indicating that the gas self-gravity becomes important 
for such models (Table [3}. Figure ^plots the initial density 
profiles for all our models. 

We consider a halo to have collapsed if within time tu, com- 
parable to the local Hubble time, the gas becomes Jeans un- 
stable within the virial extent of the DM halo. For our models 
with z = 20, 15, and 10 we set f H = 200, 300, and 500 Myr, 
respectively. We define gas to be Jeans unstable if for some 
radius r < r v ; r the virial ratio 

v = -2K(r)/W(r) (23) 

becomes less than unity. Here K(r) and W(r) are the thermal 
and potential energy, respectively, for the enclosed gas. We 
derive W(r) from the kernel-smoothed SPH gas density g: 



W(r) = -4irG [ gM(r)rdr (24) 
Jo 

iBinnev & Tremainell 19871 p. 68). For consistency, the en- 
closed gas mass M(r) in the above equation is also derived 
from the smoothed gas density g: M(r) = 4ir J gr 2 dr. In equa- 
tion J24> we ignore the potential energy of the DM halo as 
it becomes significantly smaller than the potential energy of 
the gas at advanced stages of the gas collapse. Using the 
smoothed gas density g results in W(r) being still smaller. As 
a consequence, our Jeans collapse criterion is a conservative 
one, as in a real system the potential energy for the enclosed 
gas would be larger, resulting in smaller virial ratio v. 

3. RESULTS OF SIMULATIONS 

When the flux of the external Lyman- Werner radiation is 
equal to zero, virtually all of our models collapse within one 
local Hubble time. The only exception is the least massive 
model with m vir = 10 6 M Q (r vir = 800 K) and 8=\. For the 
rest of the models, the collapse time f co u ranges from 0.02?h 
to ~ 0.9fe (see Table 0. The radius of the region becoming 
Jeans unstable, r v i r , is always slightly larger than the soften- 
ing length e, indicating that the collapse is severely under- 
resolved in our simulations. We will address the issue of nu- 
merical convergence at the end of this section. The mass of 
the unstable gas is between 3 x 10 3 and 6 x 10 4 M Q initially, 
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but becomes much larger (more than 10 6 M Q in many cases) 
by the end of the simulations due to ongoing accretion of the 
gas on the central core. In many collapsed models we could 
not continue simulations till the very end as these runs were 
virtually brought to a stop due to the extremely large gas den- 
sities developed at the center of the halo. As a result, the m V ir,2 
values we list in Table[3]should be considered as a lower limit. 
The precise value is probably not very important, as in real 
systems a star burst should occur some time after the core be- 
comes Jeans unstable, which would dramatically change the 
following evolution of the mini-galaxy. 

The main results of our work are summarized in Figure |2] 
and Table In Fi gure |21 we plot the critical value of the 
Lyman- Werner flux F2x lC tit for all our models as a function 
of m v ir, 5, and z. Models with Fn < ^2i,crit collapse within 
a local Hubble time, and do not collapse for larger fluxes. 
From Figure|2]one can see that the impact of external photo- 
dissociating radiation is most significant for halos with m w i T < 
3 x 10 6 M Q and S < 10, where even a small flux F 2 x ~ 0.01 
can prevent the halo from collapsing at a redshift of 15. As 
Figure |2j) shows, .F2i,crit is a very strong function of z. At 
z = 20, the halo with m vlI = 10 7 M© will collapse even outside 
a cosmological filament (5 = 1), whereas at z = 10 the halo 
with m v i r = 10 7 M Q can collapse only for <5 > 20 (if Fix ~ 1). 
Halos with a given virial mass are more immune to the exter- 
nal Lyman-Werner radiation at larger z due to a combination 
of several factors: (a) halos with the same m v j r have larger 
virial temperature at larger z (see Table|2j, (b) the average gas 
density in the universe grows as ~ z 3 (eq. 1221 ). and (c) the 
temperature of the cosmic microwave background grows as 
z+ 1, resulting in much stronger Compton cooling at larger 
redshifts (see eq. 1181 ). 

On a more detailed level, the impact of varying the Lyman- 
Werner flux on gas dynamics can be seen in Figure [3] Here 
we plot the central gas density and temperature as a function 
of time for our fiducial halo with m vlI = 10 7 M Q , z = 15, and 
d = 3 for different values of Fix- From this plot one can see 
that a larger photo-dissociation flux results in the central gas 
density growing more slowly and the central gas temperature 
decreasing more slowly. Also, the final temperature of the col- 
lapsed core increases with increasing flux, reaching ~ 800 K 
for logFix = -0.5. It takes longer for the core to become Jeans 
unstable (circles in Figure for larger F%x- When the flux 
reaches log/^i = 0, no collapse occurs within the evolution 
time of 300 Myr, though the central temperature decreases to 
~ 2000 K and the central density grows to ~ 30 cm" 3 . 

Figure |4] can be used to estimate the importance of differ- 
ent physical processes on the dynamics of the collapsing halo. 
Here we show the evolution of the lowest value of the virial 
ratio v within the virial extent of the halo for our fiducial case 
with m vir = 10 7 M Q , z = 15, § = 3, and F 2 x =0.1. We show 
the results for a full run (solid line) and for a few additional 
simulations with some of the physical processes turned off. 
The first most obvious effect is that turning the H2 cooling off 
changes the dynamics dramatically: when the H2 cooling is 
included, the core collapses after ~ 140 Myr, whereas when 
Ah 2 = 0, the halo is very far from a collapse even at the end of 
the simulations (after 300 Myr). Another interesting result is 
that the impact of the photo-detachment of H" (reaction 23 in 
Table[0 eq. J4]) on the overall collapse dynamics is negligible 
- despite the fact that we use the most conservative power law 
exponent a = -1 to describe the spectrum of the external radi- 
ation (corresponding to radiation produced by quasars). In a 



more realistic case, at least a fraction of the external radiation 
flux should come from stellar sources, further reducing the 
importance of the H~ photo-detachment reaction on dynam- 
ics of mini-galaxies in defunct cosmological H II regions. It 
follows that the only important photo-reaction for our objects 
is the direct photo-dissociation of H2 molecules by Lyman- 
Werner photons. 

From the analysis of Figure|4]it appears that the next most 
significant process (after H2 cooling) is the atomic hydrogen 
line cooling (long-dashed lines), followed by Compton cool- 
ing and hydrogen recombination cooling. The rest of the cool- 
ing processes have negligible impact on the evolution of the 
halo. We draw reader's attention to the fact that the impor- 
tance of the H line cooling is a strong function of the assumed 
initial temperature of the fully photo-ionized gas. E.g., low- 
ering To from 10 4 to 8 x 10 3 K reduces the initial values of 
Ahi by a factor of ~ 20 (see eq. 1141 ). It also lowers the ra- 
tio 7o/r v i r , so the overall dynamics of the collapse should not 
be affected much. We conclude that at the redshift of ~ 15 
atomic hydrogen line cooling, Compton cooling, and hydro- 
gen recombination cooling are of roughly comparable impor- 
tance for the dynamics of our mini-galaxies. 

In Figure |5] we plot the radial profiles for our fiducial halo 
for three different times t ~ 6, 144 (when the core becomes 
Jeans unstable), and 300 Myr. As one can see in Figure |5^, 
initially the core becomes unstable at the very center of the 
halo, later involving the whole virial extent of the halo. Fig- 
ure |5ji demonstrates that at the late stages of the collapse a 
dominant fraction of the total gas mass inside the virial ra- 
dius of the halo has been accreted by the core, and that the 
accretion proceeds with increasingly larger gas radial infall 
speeds. A very interesting effect can be seen in Figures |5J; 
and |5jl: after the core becomes Jeans unstable, the infalling 
gas becomes increasingly hotter near the core due to adiabatic 
compression, leading to an increasingly larger hydrogen ion- 
ization fraction near the core, and, as a result, to increased H2 
production in this area. This, in turn, increases cooling of the 
gas in the central parts of the halo, accelerating the overall col- 
lapse process. It thus appears that gas in mini-galaxies form- 
ing in the defunct cosmological H II regions can be effectively 
cooled by H2 molecules formed due to high non-equilibrium 
abundance of free electrons not only initially, before the on- 
set of the collapse, but also at the advanced stages of the col- 
lapse. This leads to a much more significant positive feedback 
of ionizing radiation on star formation in small galaxies in the 
early universe. 

From the analysis of Figure|5]it is obvious that the collapse 
is severely under-resolved in our models. E.g., one can see 
the gas infall velocity discontinuity around radius of ~ 1 3e 
(Figure^), and dramatic changes in temperature and density 
slopes at the same radius (Figure|5J;). The important question 
is then how robust are our above conclusions given that the 
collapse is not fully resolved? 

Unfortunately, it is not feasible to run higher resolution sim- 
ulations with our scalar version of the code Hydra. Instead, 
we address the above issue by re-running our models with 
m v i r = 10 7 M Q and z = 15 with two different number of parti- 
cles, 2 x 32 3 ~ 66,000 and 2 x 64 3 ~ 520,000, with the size 
of the periodic computational box being twice smaller than 
in the original simulations: 2.5 instead of 5 kpc. In half-box 
simulations the initial hydrostatic gas distribution and over- 
all collaspe dynamics is somewhat different from the full-box 
simulations, so these new runs are not meant to be directly 
compared with the original simulations. Instead, one has to 
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FIG. 2. — Critical Lyman- Werner flux i^l.crit for our models. Areas above the curves correspond to halos which do not collapse within one local Hubble time, 
and below the curves to those which do collapse, (a) Halos at fixed redshift z = 15. (b) Halos with fixed virial mass m v j, = 10 7 Mq. 
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FIG. 3. — Evolution of the gas properties at the center of the halo with 
m vir = 10 7 Mq, z = 15, <5 = 3, and logF 2 l = -co, -2, -1,-0.5,0. We show the 
gas temperature T (solid lines) and density n (dashed lines). Numbers next 
to the curves mark the values of logF9l- Filled circles (for solid lines) and 
empty circles (for dashed lines) mark the time when the halo core becomes 
Jeans unstable. 
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FIG. 4. — Evolution of the lowest value of the virial ratio v for a halo 
with m v i, = 10 7 Mq, z = 15, 5 = 3, and F21 = 0.1. Solid line: full chemistry. 
Dotted line: model with no photo-detachment of H~ (^3 = 0). Short-dashed 
line: model with Ac = (the case with A rec = is virtually identical). Long- 
dashed line: model with Ahi = 0. Dash-dotted line: model with negligible 
H2 cooling (F2i=l). 



compare the half -box low resolution simulations (which have 
the same spatial resolution as our original runs: e = 8.7 pc) 
with the half-box high resolutions simulations (spatial resolu- 
tion is twice better: e = 4.3 pc). 

The details of the half-box simulations are provided in Ta- 
ble|4](analog of the Table[3}- One can see from this table that 
in the higher resolution simulations the collapse takes place 
~ 25% earlier at a twice smaller radius, and involves initially 
2-3 times less gas. The dynamics of the collapse for both 
cases can be followed in Figure [S] where we show the evolu- 
tion of the central gas density and temperature for our fiducial 
halo. As one can see in Figure[6] the central gas density grows 



faster and reaches larger values in the higher resolution sim- 
ulations. Importantly, in the higher N case the growth rate 
of the gas density is not affected by spatial resolution even 
for some time after the onset of the Jeans instability (empty 
circles in Figure |6j, whereas in the lower resolution case the 
sharp change in the density growth rate takes place at the on- 
set of the instability. This indicates that at least the initial 
stage of the core collapse is resolved in our higher resolution 
simulations. Interestingly, the final central temperature of the 
gas is the same in the both cases (see Figure[5}. 

It is easy to see why in the higher resolution case the col- 
lapse takes place earlier. Indeed, from Figure [2 one can see 
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FIG. 5. — Radial profiles for a halo with m v j r = 10 7 Mq, z= 15, S = 3, and logKi = - 1 at three different times: at the early epoch (1 = 6 Myr, dash-dotted 
curves), around the time when the gas became Jeans unstable (r = 144 Myr, solid curves), and at the end of the simulations (; = 300 Myr, dashed curves). Vertical 
dash-dotted lines correspond to the spatial resolution 1.3e. (a) Virial ratio v. Parts of the curves below the horizontal line logf = are Jeans unstable, (b) Mass 
M of the gas enclosed within radius r (thick lines) and radial gas velocity V,- (thin lines). Horizontal solid line corresponds to V r = 0. (c) Gas density n (thick 
lines) and temperature T (thin lines), (d) Fractional abundances / H + (thick lines) and /h 2 (thin lines). 



TABLE 4 

Half-Box Simulations for a Halo with logm vir = 7 and z = 15 



N 


.5 


logF 2 i 


n c 


rtc/nout 


logftl.crit 


fcoll 


'•vir.l 




F vir,2 


»M v jr,2 




io~ 21 


ergs s -1 cirT 2 sr~' Hz~ 


cm 3 






Myr 


pc 


M 


pc 


M Q 


2 X 32 3 


1 


-00,-3,-2,-1.5-1 


5.8 X 10~ 3 


6.9 


-1.25 


103 


10 


9.0 x 10 3 


10 


2.4 X 10 5 


2 X 32 3 


3 


-00,-2,-1,-0.5,0 


0.018 


6.9 


-0.75 


51 


11 


1.0 x 10 4 


10 


7.5 x 10 5 


2 X 32 3 


10 


-00,-1,0,0.5,1,00 


0.062 


7.2 


0.25 


31 


11 


1.2 x 10 4 


10 


2.4 x 10 6 


2 X 32 3 


30 


-oo,-l,0,l,oo 


0.22 


8.0 


CO 


21 


26 


5.0 x 10 4 


10 


6.1 x 10 6 


2x64 3 


1 


-00,-3,-2,-1.5,-1 


7.8 x IO -3 


9.2 


-1.25 


85 


5 


4.4 x 10 3 


5 


5.7 x 10 4 


2x64 3 


3 


-00,-2,-1,-0.5,0 


0.024 


9.2 


-0.75 


40 


5 


4.9 x 10 3 


5 


1.4 x 10 5 


2x64 3 


10 


-00,-1,0,0.5,1,00 


0.083 


9.6 


0.25 


22 


6 


6.3 x 10 3 


5 


3.4 x 10 5 


2 x 64 3 


30 


-oo,-l,0,l,oo 


0.30 


11 


oo 


15 


9 


1.4 x 10 4 


5 


9.0 x 10 5 



NOTE. — Here N is the total number of particles (DM + gas). The rest of the notations are the same as in Tablel3l 
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FIG. 6. — Evolution of the gas temperature (solid lines) and density (dashed 
lines) at the center of the halo with m v j r = 10 7 Mq,z = 15, 8 = 3, and logf^l = 
-1 in half-box simulations. Thick and thin lines correspond to high (N = 2 X 
64 3 ) and low (N = 2 x 32 3 ) resolution simulations. Filled (for solid lines) and 
empty (for dashed lines) circles mark the time when the halo core becomes 
Jeans unstable. 



that all our initial hydrostatic gas density profiles have a cusp 
with a small slope at the center of the halo. As a result, the 
higher the spatial resolution of the simulations, the larger is 
the initial density of the resolved gas in the core, leading to 
stronger cooling and faster evolution. As the slope of the 
density cusp is small, we expect this process to converge at 
some large N. Indeed, as we progress to larger N, a very 
moderate increase in the resolved gas density (and hence of 
the cooling rate) would correspond to much smaller mass of 
the gas affected by the cooling, resulting in negligible impact 
on the dynamics of the system for large enough N. More 
quantitatively, if the initial gas density near the halo cen- 
ter is oc r~ a and the spatial resolution is e, then the mass 
of the unresolved core is M c oc e 3 " Q /(3-aO and the average 
core density is p c oc e~ a /(3 — a). As the cooling per atom 
is proportional to p c , the global core cooling is proportional 
to M c p c oc e 3_2a /(3 - a) 2 , which approaches zero as e — ► 
(N — > oo) for shallow cusps with a < 1.5. All of our models 
initially had a = 0.11 . . .0.55. 

Figure allows a detailed comparison of the final gas pro- 
files in our half-box simulations for the cases of low and high 
resolution. From Figure 05 one can see that the final mass 
of the collapsed gas is comparable in the both cases, and that 
the gas infall velocity discontinuity is even sharper and takes 
place at a twice smaller radius in the higher N case. Most 
importantly, the higher resolution simulations confirm and re- 
enforce our above conclusion that the non-equilibrium forma- 
tion of H2 molecules in our objects is not restricted to the ini- 
tial stage of the collapse and is instead an ongoing process, 
becoming potentially even more important at the advanced 
stages of the collapse. This is evident in Figures Q; and 01, 
which show that in the higher resolution simulations the gas 
temperature near the core is even higher, with the hot gas zone 
extending to even smaller radius, than in the lower N case, 
resulting in higher non-equilibrium electron abundance and 
hence higher H2 abundance in the core region. 

It is not clear from the above results if our higher resolu- 



tion half-box simulations are close to numerical convergence. 
Nevertheless, by comparing our low and high N simulations 
we can make a robust conclusion that the models presented 
in this paper capture the essence of the gas collapse process 
in primordial mini-galaxies forming in defunct cosmological 
H II regions. Our results are conservative in the sense that 
higher resolution, and presumably more realistic, simulations 
would lead to mini-galaxies forming even faster and poten- 
tially even under larger fluxes of external Lyman-Werner ra- 
diation. 

4. DISCUSSION AND CONCLUSIONS 

The traditional point of view is that soon after the for- 
mation of the first stars in the universe inside small ha- 
los at z > 20 the formation of the next generation of mini- 
galaxies with T v i r < 10 4 K was suppressed due to build 
up of the m etagalactic H2 photo-dissociating background 
jHaiman et alJ |2000|). Recent cosmologic al ( Ricot tiet alJ 
120021) and semi-cosmological (Q'She a~et al J 120051) simula- 
tions with radiative transfer questioned the above argument 
by showing that in defunct cosmological H II regions in over- 
dense regions of the early universe (z > 10) the formation 
of the mini-galaxies can be boosted due to non-equilibrium 
chemistry. In our paper we study the above mechanism of 
positive feedback of ionizing radiation on small galaxy for- 
mation by running a large set of simulations of galaxy for- 
mation inside defunct H II regions. By employing a non- 
cosmological approach we produced a sequence of models 
with precisely controlled and reproducible (though idealized) 
initial and boundary conditions. This helped us to focus on 
the most relevant physical processes at work in our objects. 

Most importantly, we derived the critical fluxes of the back- 
ground Lyman-Werner radiation sufficient to prevent the col- 
lapse of gas in mini-galaxies with T W1I < 10 4 K at z = 10-20. 
We showed that the positive feedback mechanism becomes 
much more efficient at larger z fo r galaxies with give n virial 
mass, and confirmed the result of Ricotti et al. (2002) that it 
is efficient inside filaments with the overdensity of > 10. We 
studied the details of the process of the gas collapse in a pre- 
ionized mini-galaxy exposed to external Lyman-Werner radi- 
ation. We conclude that, similarly to the canonical case of 
small galaxy formation in non-ionized regions of the universe 
( Mach acek et alJl2001h . the reaction of photo-detachment of 
H" has a negligible impact on the dynamics of the collapse. 
This leaves the direct photo-dissociation of H2 molecules by 
metagalactic Lyman-Werner radiation as the only important 
photo-reaction. We identified the H2, hydrogen line, Comp- 
ton, and hydrogen recombination coolings as the only impor- 
tant cooling processes. Hydrogen collisional ionization cool- 
ing and Bremsstrahlung cooling are shown to be not important 
for our objects. 

We found that for the mini-galaxies forming in defunct H II 
regions the non-equilibrium H2 formation has an important 
dynamical effect not only during the initial phase leading to 
the core of the galaxy becoming Jeans unstable, but also later 
on during the advanced stages of the core collapse. This 
mechanism allows the core of the mini-galaxy to accrete sig- 
nificant mass by the time the central star burst disrupts the 
accretion. 

Our results give further support to the iRicotti et alJ |2002) 
picture of formation and evolution of mini-galaxies in the 
early universe. In this picture, first mini-galaxies form at high 
redshift (z > 20), when the metagalactic Lyman-Werner back- 
ground is close to zero, due to cosmological non-zero resid- 
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FIG. 7. — Final radial profiles for a halo with m v j r = 10 7 Mq, z = 15, 8 = 3, and logF2i = _ 1 for half-box simulations at two different resolutions: N = 2 X 64 3 
(solid lines) and N = 2 X 32 3 (dashed lines). Vertical lines mark the corresponding spatial resolutions, 1.3e. (a) Virial ratio v. Parts of the curves below the 
horizontal line \ogu = are Jeans unstable, (b) Enclosed gas mass M (thick lines) and radial gas velocity V,- (thin lines). Horizontal solid line corresponds to 
V r = 0. (c) Gas density n (thick lines) and temperature T (thin lines), (d) Fractional abundances / H + (thick lines) and /h 2 (thin lines). 



ual fraction of free electrons resulting in efficient H2 cooling. 
The stars in the mini-galaxy ionize gas inside the halo and in 
nearby halos, creating conditions for efficient H2 cooling due 
to non-equilibrium electron abundance after the Pop III stars 
(believed to be very massive and short-lived) die. As we con- 
firmed here, this process is tolerant to relatively large fluxes 
of the external Lyman-Werner radiation. The process can re- 
peat itself with a time scale of a few t ens of millions years 
(iRicotti et alJl2002t lO'Shea et alJl2005t this paper), allowing 
the mini-galaxy to survive through the "difficult times" of the 
universe being filled with the H2 photo-dissociating radiation 
until the universe is reionized before z ~ 7. These repetitive 
local star bursts will result in the metallicity of the local in- 
tergalactic medium increasing with time. If the metallicity 
becomes larger than ~ 10~ 3 Zq, the star-formation will switch 
from Pop III to Pop II mode (e.g. Bromm et al. 2001p, pro- 
ducing a long-living stellar population, which should survive 



until the present time. The size of the collapsed cores in our 
simulations is very small (< 5 pc) and is comparable to the 
sizes of globular clusters. It is tempting then to associate this 
final Pop II star burst at the core of the mini-galaxy with the 
formation of a metal-poor globular cluster - at least for the ha- 
los with the l owest values of the spe cific angular momentum. 
We showed dMashchenko & Sillsll2005bl) that the DM halos 
of such objects can be tidally stripped when they are accreted 
by larger galaxies. This could be considered as a revised ver- 
sion of the primo r dial pi cture of globular cluster formation of 
iPeebles & Dickd i 19681) and iPeeblesI (1984). Alternatively, 
these surviving mini-galaxies could end up contrib uting to 
the po pulation of dwarf spheroidal galaxies jRicotti & Gnedinl 
2005). 

In § we discussed the impact of the limited resolution 
on the results of our simulations. Our conclusion was that 
even though our highest spatial resolution simulations could 
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not completely resolve the dynamics of the collapsing core, 
we believe that we captured the essence of the collapse pro- 
cess. We want to emphasize that in our simulations the 
physics driving the collapse (primordial chemistry and cool- 
ing) is resolved much better than hydrodynamics. As Fig- 
ure |5Jl demonstrates, during the most critical initial phase of 
the core collapse the enhanced non-equilibrium production of 
H2 molecules is not restricted to the core region, and instead 
takes place over a wide range of distances from the center of 
the halo. Another confirmation of the chemistry being well re- 
solved in our models can be found in Table|4] where one can 
see that the critical fluxes of the external photo-dissociating 
radiation F2i.ciit are the same for both lower and higher reso- 
lution (half- box) simulations. 

Recently, Glover, Savin, & Jappsen (2005) used the results 
of their simulations of the collapse of gas in a photo-ionized 
small galaxy in the early universe to argue that the uncertain- 
ties in the rate coefficients for the reactions 8 and 16 (see Ta- 
ble Q can lead to an order of magnitude uncertainty in the 
value of the critical flux of the external photo-dissociating 
radiation. We cannot comment on this as we have not ex- 
plored in our models the consequences of the uncertainties in 
the rates for the ab ove reactions being large. If the results of 
iGlover et al J (12005) are confirmed, this would potentially lead 
to up to a factor of a few error in our values of /^l.cnt, with all 
the values being either systematically larger or smaller. We 
want to emphasize that it would not change our qualitative re- 
sult that the efficient formation of mini-galaxies in fossil cos- 
mological H II regions is possible even under relatively large 
levels of the metagalactic Lyman- Werner radiation. 

How our results are affected by the simplifying assump- 
tions we have made? One of our model assumptions is that 
the process of photo-evaporation of the mini-galaxy is com- 
plete, with the gas becoming 100% ionized and reaching the 
state of hydrostatic equilibrium in the halo gravitational po- 
tential. This assumption allowed us to significantly reduce 
the uncertainty in the initial conditions. In reality, in many 
(perhaps most) cases, the gas in our objects would only be 
partially photo-ionized. The reason for this is that the typi- 
cal sources of ionizing radiation at this early epoch would be 
massive stars (individual or in small clusters), with their life- 
times being too short to allow the gas in the galaxy, which is 
being photo-evaporated, to reach a state of hydrostatic equi- 
librium. Assuming that initially the gas was almost neutral 
and had a temperature commensurate with the virial tempera- 
ture of the halo T vlr < 10 4 K, an incomplete photo-evaporation 
would result in a partially ionized gas with the core density 
larger than in the fully photo-ionized state. It is reasonable to 
assume that such a system should collapse at least as easily 
(or even easier) as a fully photo-ionized one, as long as the 
electron abundance of the photo-ionized gas is significantly 



larger than the residual cosmological value of ~ 2 x 10 . In 
this sense, the above assumption is a conservative one. 

Our assumption of the gas being optically thin to radia- 
tion is also a conservative one, as the self-shielding of the gas 
would reduce the Lyman- Werner radiation flux near the center 
of the halo, where the production of H2 molecules takes place 
during the advanced stages of the core collapse. In our fiducial 
halo simulated with the highest spatial resolution (half-box 
simulations with N = 2 x 64 3 ), the optical depth for Lyman- 
Werner photons traveling from the outside of the halo to the 
smallest resolved radius of ~ 5 pc is around unity, indicating 
that self-shielding can be important for our objects (especially 
at the late stages of the c ore co llapse). 

Nagakura & Omukai (2005) used one-dimensional hydro- 
dynamic simulations incorporating non-equilibrium primor- 
dial gas chemistry to show that the gas cooling due to HD 
molecule can be an important dynamical factor at the late 
stages of the gas collaps e in prev i ously photo-ionized mini- 
galaxies. Recently, Lipovka et al. (2005) published a revised 
HD cooling function, with the cooling efficiency being larger 
than previously reported (especially for hot and dense gas). 
This could potentially make the HD molecule cooling even 
more important for the dynamics of the collapsing gas in small 
primordial galaxies. We did not include HD chemistry in our 
model as these reactions would only become important when 
the collapsing core is no longer well resolved in our simula- 
tions. If included, the additional cooling due to HD molecules 
would only strengthen our conclusions. 

If gas in photo-evaporated mini-galaxies has significant 
global angular momentum this could affect our results and 
conclusions. The effect of angular momentum will be 
strongest during the advanced stages of the gas collapse. In- 
stead of a spherical core, a central disk will be formed, reduc- 
ing the central gas density. We argue that the non-zero angular 
momentum should not change qualitatively our main conclu- 
sions, as the non-equilibrium H2 formation, leading to the gas 
in the mini-galaxy becoming Jeans unstable, takes place over 
a wide range of radial distances (see thin s olid line in Fig- 
ure|5Jl), not just in the central densest region. Machacek et al. 
(2001) noted that the ability of a cloud to collapse does not 
seem to strongly depend on the angular momentum of its DM 
halo. In any event, our results should be directly applicable to 
mini-galaxies with low values of the global gas angular mo- 
mentum. 



The simulations reported in this paper were carried out on 
McKenzie cluster at the Canadian Institute for Theoretical As- 
trophysics. HMPC acknowledges support from the Canadian 
Institute for Advanced Research and NSERC. 
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